Effect of polarized gluon distribution on spin asymmetries 
for neutral and charged pion production 
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A longitudinal double spin asymmetry for vr" production has been measured by the PHENIX 
collaboration. The asymmetry is sensitive to the polarized gluon distribution and is indicated to be 
positive by theoretical predictions. We study a correlation between behavior of the asymmetry and 
polarized gluon distribution in neutral and charged pion production at RHIC. 
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I. INTRODUCTION 

Determination of the polarized parton distribution 
functions (PDFs) is crucial for understanding the spin 
structure of the nucleon 0. As is known well, the proton 
spin is composed of the spin and orbital angular momen- 
tum of quarks and giuons. Several parametrizations of 
the polarized PDFs have been proposed, and have suc- 
cessfully reproduced experimental data H, H, ^ IE IS- 
In particular, the amount of the proton spin carried by 
quarks is determined well by global analyses with the po- 
larized deep inelastic scattering (DIS) data. The value is 
about AS = 0.1 ^ 0.3, whereas the prediction from the 
naive quark model is AE — 1 . This surprising result leads 
to extensive study on the gluon polarization. The cur- 
rent parametrizations suggest a large positive polariza- 
tion of gluon. However, our knowledge about the polar- 
ized gluon distribution /S.g{x, Q^) is still poor, since the- 
oretical and experimental uncertainties are rather large. 
The determination of Ag{x,Q^) gives us a clue to the 
proton spin puzzle. 

The RHIC is the first high energy polarized proton- 
proton collider to measure Ag{x,Q^) 1^. We can ex- 
tract information about Ag{x, Q^) through various pro- 
cesses, e.g., prompt photon production, jet production, 
and heavy flavor production. These processes are quite 
sensitive to A.g{x,Q^), since giuons in the initial state 
associate with the cross section in leading order (LO). 

Recently, the PHENIX collaboration has reported re- 
sults for inclusive i:^ production pp — s- tt'^X ^] which 
is also likely to be sensitive to Ag{x,Q^). The double 
spin asymmetry was measured in longitudinally polar- 
ized proton-proton collisions at RHIC in the kinematical 
ranges: center-of-mass (cm.) energy a/s = 200 GeV and 
central rapidity I77I < 0.38. The data imply that the 
asymmetry might be negative at transverse momentum 



Pt — I ^ 3 GeV. The lower bound of the 7r° asym- 
metry at low pt has been considered, and a slight nega- 
tive asymmetry by modifying Ag(x, Q^) has been demon- 
strated in Ref. [3- However, there is no theoretical pre- 
dictions indicating large negative asymmetry. 

In this paper, we study the behavior of the tt*' dou- 
ble spin asymmetry correlated with Ag{x,Q^) in SeclTTl 
By using three types of Ag{x,Q^), we suggest that the 
asymmetry in large pt region is more sensitive to the 
functional form of Ag{x, Q^)- An impact of the new data 
on determination of Ag{x) is discussed in terms of un- 
certainty of the asymmetry coming from the polarized 
PDFs. Furthermore, we discuss a spin asymmetry for 
charged pion production in Sec. IIIII An asymmetry tak- 
ing the difference of cross sections for tt^ and tt^ pro- 
duction is proposed, and it is useful to discuss the sign 
of Ag{x) in the whole px region. The Summary is given 
in Sec. ifvl 



II. 



SPIN ASYMMETRY FOR NEUTRAL PION 
PRODUCTION 

A. Ambiguity of the polarized cross section 



First, we describe the longitudinal double spin asym- 
metry for tt" production. It is defined by 



_ [da^^ — da^ ]/dpT dAa/dpx 

[da^+ + da-{ ]/dpT da/dpr 



(1) 
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where px is the transverse momentum of produced pion. 
dajih' denotes the spin-dependent cross section with def- 
inite helicity h and h' for incident protons. 

The cross sections can be separated short distance 
parts from long distance parts by the QCD factorization 
theorem. The short distance parts represent interaction 
amplitudes of hard partons, and are calculable in the 
framework of perturbative QCD (pQCD). On the other 
hand, the long distance parts such as PDFs and fragmen- 
tation functions should be determined by experimental 
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data. The polarized cross section Act is written by 



dpi 



/ dij f dxa I dxb 
xA/,(x„Q2)A/(,(xf„Q2) 



d{pT,v) 



dt 



(2) 



where Afi{x, Q^) is the polarized PDFs, and D^'\z, Q^) 
is the spin-independent fragmentation function decaying 
into pion c — > tt" with a momentum fraction z. The sum 
is over the partonic processes a + b ^ c + X associated 
with vr" production. is the Jacobian which transforms 
kinematical variables from t and z into px and 77 of the 
produced tt*^. Act describes the polarized cross section 
of subprocesses. The partonic Mandelstam variables s 
and t are defined by s = {pa + PbY and t — (pa — Pc)"^ 
with partonic momentum pi, respectively. The squared 
cm. energy s is related to s through s = XaXbS and set 
as ^/s — 200 GeV. The pseudo-rapidity 77 is limited as 
I77I < 0.38 in the PHENIX acceptance. 

In this analysis, the cross sections and the spin asym- 
metry are calculated in LO level. Rigorous analysis of 
Oipt^g) next-to-leading order (NLO) calculation has been 
established in Ref. We believe that the qualitative 

behavior of the asymmetry does not change, even if NLO 
corrections are included in our study. In numerical calcu- 
lations, we adopt the AAC set |5| as the polarized PDFs 
and the Kretzer set |0| ^ the fragmentation functions. 
We choose the scale — Pr- 

The partonic subprocesses in LO are composed of 
0{a1) 2^2 tree-level channels listed as gg q{g)X, 
19 ^ q{9)X^ m m qig, q')X, qq' -> qX, and 

qq' ~^ qX including channels of the permutation q ^ q. 
Main contribution to the polarized cross section comes 
from gg — + q{g)X and qg — > q{g)X channels with con- 
ventional PDFs and fragmentation functions. The gg 
contribution dominates in low pT region and steeply de- 
creases with pt increases. Then, the qg process becomes 
dominant in larger pT region. The crossing point of these 
contributions however depends on parametrization of the 
polarized PDFs. In both cases, the spin asymmetry for 
TT^ production is sensitive to the gluon polarization. 

As mentioned above, the partonic cross section Act is 
well-defined in the pQCD framework. Hence, as a cause 
of inconsistency with the PHENIX data, we consider the 
ambiguity of long distance parts: fragmentation func- 
tions and PDFs. 

The fragmentation into tt*^ includes all channels 
q, g, 5 — !■ 7r°. Each component of the fragmentation func- 
tions can be determined by global analyses with 
several experiments |0 . The unpolarized cross sec- 
tion measured by the PHENIX are consistent with 
NLO pQCD calculations within model dependence of 
. These precise measurements give strong constraint 



on the fragmentation functions. Significant modification 
of them would not be expected. Therefore, the fragmen- 
tation functions are not the source of the negative asym- 
metry even if they have uncertainty to some extent. 

In the polarized reaction, kinematical ranges and the 
fragmentation functions are the same as the unpolarized 
case except the polarized PDFs. For the polarized quark 
distributions Aq{x) and Ag(a::), the antiquark distribu- 
tions and their flavor structure are not well known. For 
TT*^ production, subprocesses are (light quark) flavor blind 
reaction, and the predominant qg process depends on the 
sum /S.q{x) + Ag(a;) which is relatively determined well 
by the polarized DIS data 0, and so ambiguities of the 
polarized quark distributions can be neglected. In con- 
sequence, the undetermined polarized gluon distribution 
Ag{x) remains as the source of the uncertainty of the 
asymmetry. 



B. Correlation between the spin asymmetry and 
the polarized gluon distribution l^g{x) 

To investigate a role of Ag{x) for the behavior of the 
asymmetry, we prepare three functional forms as shown 
in Fig^ Solid curve shows Ao(a;) by the global analysis 
with the polarized DIS data ^ . Dashed and dot-dashed 
curves show two artificial modified A(7(x), respectively. 
The sample-1 distribution has a node. The gluon dis- 
tribution with a node has been indicated in the paper 
by Jager et ai, Q. Our distribution is negative in the 
small X region, and positive in the large x region. It 
has opposite signs of Ag{x) shown in Fig. 2 of their 
paper. The sample-2 distribution is small negative in 
the whole x region. Their distribution is similar to the 
sample-2 rather than the sample-1. It shows barely pos- 
itive at small x, while the sample-2 is negative. Since 
the sample-1 and 2 are within the Ag{x) uncertainty by 
the AAC analysis, these distributions can be adopted as 
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FIG. 1: Polarized gluon distributions Ag{x) at pr = 2.5 
GeV. Solid, dashed, and dot-dashed curves indicate the AAC, 
sample-1, and 2 distributions, respectively. 
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FIG. 2: Spin asymmetries for tt" production by using three 
different Ag{x) in Fig. 



a model of Ag{x). These are taken account of the 
dependence by the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equation with the polarized quark and 
antiquark distributions. 

We discuss the behavior of the spin asymmetry asso- 
ciated with the functional form of Ag{x). The obtained 
asymmetries with these giuon distributions are shown in 
Fig. 121 We find that the asymmetry for the AAC Ag{x) 
is positive in the whole px region. The asymmetries for 
the sample-1 and 2 become negative at low pT- In partic- 
ular, we obtained the negative asymmetry in the whole 
Pt region by using the sample-2 Ag{x). Furthermore, 
one can see variations of these asymmetries at large pt- 

The asymmetry for the AAC is positive and increases 
with Pt- The positive polarization for Ag{x) generates 
positive contributions of gg and qg processes which dom- 
inate in the 7r° production. In this case, the asymmetry 
cannot become negative. 

The positive Ag{x) is suggested by the recent global 
analyses with the polarized DIS data 0, |3, 0, 0, U- Al- 
though these analyses obtain good agreement with the 
experimental data, the Ag{x) cannot be determined and 
it has large uncertainty. Therefore, we cannot rule out 
the negative polarization for Ag(x). There is a possibility 
of the negative asymmetry with the modified Ag{x). 

For the sample-1 in Fig. |21 the asymmetry is slight 
negative at low px and changes into positive at px = 5 
GeV. As is mentioned in Ref. 0], the Ag{x) with a node 
has the possibility of making the small negative asym- 
metry at low Pt- In the region px < 5 GeV, we find 
that contributions of gg and qg processes are negative, 
respectively. To make negative gg contribution would 
be needed opposite polarizations of Ag{x) at Xa and Xb- 
Computed by using several shapes of Ag{x) with a node, 
the gg contribution is not always negative. The contri- 
bution basically depends on the shape of Ag{x) even if 
it has a node. 

In the region > 3 GeV, the gg contribution changes 
into positive, and dominates in the region px < 10 GeV. 



This is because that the node rapidly shifts toward low- a; 
direction due to evolution with px- Therefore, the 
positive polarization for Ag(x) at medium x contributes 
predominantly to the positive asymmetry via the gg pro- 
cess. Furthermore, the asymmetry at large is sensitive 
to the behavior of Ag(x) at medium x. 

As another possibility of the negative asymmetry, we 
choose slight negative polarization for Ag{x). In this 
case, the gg contribution is positive while the qg con- 
tribution is negative. The asymmetry is determined by 
the difference between two contributions. The gg and 
qg contributions are proportional to (Ag)^ and Ag, re- 
spectively. The gg contribution is more sensitive to the 
behavior of Ag(x). In particular, the behavior at low x 
significantly affects on the contribution at low px since 
the value of x""" in Eq. Q is rather small. In order 
to make the positive gg contribution smaller, the Ag{x) 
for the sample-2 is taken small polarization at low x as 
shown in Fig. ^ 

In Fig. 12 as far as the sample-2 is concerned, the asym- 
metry indeed becomes negative in the whole px region. 
In the region px < 3 GeV, the small negative polar- 
ization for Ag(x) generates slight positive contribution 
of the gg process. In this case, the gg contribution is 
the same order of magnitude as the qg contribution, and 
almost cancel out the negative contribution. The asym- 
metry is therefore determined by other processes. The 
total contribution of the processes except above two pro- 
cesses becomes slight negative. Above the region, the gg 
contribution rapidly decreases with px increases. The qg 
process becomes dominant contribution, which provides 
the negative asymmetry Thus, the negative asymme- 
try can be obtained in the whole px region by using the 
negative Ag{x) which makes the qg contribution larger 
than the gg contribution. 

In the sample-2, we should note that the magnitude 
of Ag(x) at the minimum point cannot be large. This 
is because that the shape of Ag{x) is rapidly varied by 
the evolution, the minimum point of Ag(x) shifts to- 
ward low- a; and the width broadens. At moderate px, 
the gg process is more sensitive to the low- a; behavior of 
the evolved Ag{x) than the qg process. If the Ag{x) is 
taken large negative polarization at the minimum point, 
the magnitude of the gg contribution becomes rapidly 
large compared with the qg contribution, and then the 
asymmetry becomes positive at moderate px- The small 
negative Ag(x) is therefore required to obtain the nega- 
tive asymmetry in the whole px region. 

In above two cases at low px, we cannot also obtain 
negative value exceeded the lower limit —0.1% that is 
suggested in Ref. 0. Furthermore, even if the asym- 
metry is positive, the magnitude is below 1%. As we 
discussed, the functional form of Ag(x) needs some re- 
straints to make the asymmetry negative. It is difficult 
to obtain sizable value in comparison with the positive 
case. 

At large Pt, the difference of the obtained asymmetries 
remarkably reflects the medium-x behavior of Ag{x). Ex- 
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perimental data in the region is useful to determine the 
^g{x). For instance, the asymmetry for the sample-2 be- 
comes rather larger to negative direction. If future pre- 
cise data indicate the negative asymmetry in the region, 
the ^g{x) requires significant modification of its func- 
tional form. It has the potential of the negative gluon 
contribution to the nucleon spin. In order to understand 
the behavior of ^g{x) in detail, we require experimental 
data covering a wide pr region. 



C. Uncertainty of the spin asymmetry 

Next, we consider the effect of the 7r° data on the l^g{x) 
determination in terms of the uncertainty estimation for 
the spin asymmetry. The large uncertainty of i^g{x) im- 
plies the difficulty of extracting the gluon contribution 
from the polarized DIS data. We are therefore interested 
in constraint power of the new data on lS.g{x). If the 
experimental data are included in a global analysis, the 
asymmetry uncertainty will be bounded within statisti- 
cal error range. See, for example. Fig. 2 of Ref. [2]. As 
far as evaluation of the constraint is concerned, the un- 
certainty can be compared with statistical errors of the 
data, although it is rough evaluation. 

The asymmetry uncertainty coming from the polar- 
ized PDFs is defined by a polarized cross section uncer- 
tainty divided by a unpolarized cross section: SAJ^j^ = 
SAa"^" I a'^^' . The cross section uncertainty is obtained by 
taking a root sum square of uncertainties of all subpro- 
cesses. These uncertainties are estimated by the Hessian 
method, and are given by 



i,3 



(3) 

where k is the index of subprocesses. is a optimized 
parameter in the polarized PDFs. Hij is the Hessian 
matrix which has the information of the parameter er- 
rors and the correlation between these parameters. The 
Ax^ determines a confidence level of the uncertainty, and 
is estimated so that the level corresponds to the Icr stan- 
dard error. We choose the same value as the AAC anal- 
ysis Further, the gradient terms for the subprocesses 
dAa^ ) / dai are obtained by 



dAal 
dpT 



E 

a.h.c 



drj 



dx„ 



dxh 



Afb[Xb) + A}a{Xa) 



dai 
d{i, z) 
d{pT,v) 



A&- 



ab—^cX 



is,i) 



dt 



da. 



D: (z), (4) 



The gradient terms for the polarized PDF dAfa{xa)/dai 
are analytically obtained at initial scale Qqi f^nd are nu- 
merically evolved to arbitrary scale by the DGLAP 
equation. 




FIG. 3: Comparison of the asymmetry uncertainty SAJ^j^ 
with the statistical errors for ^/s = 200 GeV. 



In Fig. 13 the asymmetry uncertainty is compared 
to the statistical errors of the experimental data by the 
PHENIX 8]. In this comparison, we exclude the data 
at pt = 1-5 GeV. This is because that the data might 
have contribution from soft physics, and it might not be 
explained as physics of a hard process. We have no idea 
whether such data can be included in the global analy- 
sis. From this figure, we find that the uncertainty almost 
corresponds to the experimental errors, and is mainly 
composed of the uncertainty of Ag{x). This fact indi- 
cates that the present tt" data have the same constraint 
on Ag{x) as the polarized DIS data. At this stage, one 
cannot expect to reduce the Ag{x) uncertainty even if 
these data are included into the global analysis. How- 
ever, the asymmetry uncertainty is very sensitive to the 
Ag(x) uncertainty. The 7r° production has the potential 
to become a good probe for Ag{x) by future precise data. 

It should be noted that symmetric uncertainty is shown 
in order to compare with the statistical errors in Fig. |31 
The lower bound is however incorrect because a lower 
limit of the asymmetry is not taken into account. As 
mentioned in previous subsection, the asymmetry can- 
not exceed —0.1 % at low pr where the gg process domi- 
nates. Although asymmetric uncertainty should be esti- 
mated, such uncertainty cannot be obtained by the Hes- 
sian method. We therefore need further investigation of 
the lower bound for the asymmetry uncertainty. 



III. SPIN ASYMMETRY FOR CHARGED PION 
PRODUCTION 

We discuss the spin asymmetry for charged pion pro- 
duction, TT"*" and TT^. Unpolarized and polarized cross 
sections can be similarly calculated by using the frag- 
mentation functions decaying into charged pion D'" in 
Eq. ((SJ. We show asymmetries with the AAC Ag{x) 
and sample-2 Ag(x) in Fig. 01 In the asymmetries for 
the AAC Ag{x), one can see differences among them in 
large px region where the qg process is dominant. The 
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FIG. 4: Asymmetries for neutral and charged pion produc- 
tions with the AAC and sample-2 Ag{x) sets. 



polarized cross sections of qg process for 7r+ and tt pro- 
duction are written by 



99 



= Ag 




(5) 



where the symbol (g> denotes convolution integral in Eq. 

z indicates the quark flavor, and is taken as i — 
u, d, s, u, d, and s. Actual calculation includes permu- 
tated terms of Xa and Xb- There are following relations 
among the fragmentation functions for charged pion: 



1 -^9 ' 



(6) 



and the fragmentation functions for neutral pion are de- 
fined by 



A" =(A" +Df)/2 



(7) 



For -K^ production, the contribution associated with 
Au distribution is enhanced by the fragmentation func- 
tion DJJ^ . Increasing asymmetry for 7r+ production is 
caused by positive contribution from Au distribution, 
whereas decreasing asymmetry for tt^ production comes 
from negative Ac? distribution. 

On the other hand, the asymmetries for the sample-2 
Ag(x) are almost the same. The differences among them 
depends on the magnitude of lS.g{x), since the asymmetry 
is proportional to A.g{x) as written in Eq. ((S)). If the 
absolute value of A.g{x) is small, there are not significant 
differences among the asymmetries for tt*^, tt"'', and tt^ 
productions. 

In order to determine A.g{x) with its sign by using 
charged pion production, let us propose an interesting 
observable which is defined by 



A 



LL 



(8) 



The behavior of the asymmetry is sensitive to the sign of 
Ag{x) because the contribution of the gg processes are 
eliminated and one of the qg process becomes dominant 
in the whole px region. The polarized cross section for 
gg ~^ gg process is given by 



Aa^g^ =Agr^Ag^Df . 



)Act9 



(9) 



This contribution is cancelled out due to — . 
For the same reason, gg —^ qq process does not also con- 
tribute by summing fragmentation functions for flavors: 
J2i Df ~ J2i ■ As the similar case, the contribu- 
tions of qq gg, qq q'q' processes are also vanished. 
The unpolarized cross section can be similarly calculated 
with unpolarized PDF's and partonic cross sections. 

The asymmetry can be obtained by the difference of 
qg process. The second term of Eq. (0) is cancelled out 
for the same reason of gg — > gg process. And then, the 
asymmetry is consequently given by 



Ag » {Auy -Ady)(s {Dl - D^)ra Aaia- 



(10) 

where Afv{~ Af — A/) is a polarized valence quark dis- 
tribution. The following relations among the fragmenta- 
tion functions are assumed by the isospin symmetry. 



= Df - Df = Df 



Dl 

^2 



(11) 



This relation is used in parametrization of the fragmen- 
tation functions jl 11] . 

In the asymmetry in Eq. ambiguity of fragmen- 

tation function is removed by the cancellation of the 
convolution part. Another ambiguity from the fragmen- 
tation functions can be also cancelled between numerator 
and denominator. In addition, Au-u — Ady is determined 
well, since its first moment is constrained by neutron and 
hyperon beta decay constants |0,l3||, iSl] ■ Of course, un- 
polarized PDF's are precisely determined in comparison 
with the polarized PDF's. This asymmetry can be de- 
fined by well known distributions without Ag; therefore, 
we can effectively extract information about Ag(x) in- 
cluding its sign. 

Figure. shows the asymmetry defined by Eq. ((SJ. 
Solid and Dotted curves are asymmetries with AAC 
Ag{x) and —Ag(x). We find large asymmetries in both 
cases. In particular, the asymmetry with —Ag{x) is neg- 
ative and the absolute value is large in comparison with 
single pion production. Since the asymmetry is domi- 
nated by qg process in the whole px region, the differ- 
ence of the sign of Ag{x) is markedly refiected in the 
asymmetry. 

We mention the contribution of qq process to the 
asymmetry. In the region 8 < < 13 GeV, the qq 
contribution accounts for 10-15% of the polarized part 
{Aa'^ — Aa-"^ ), and 27-56% of the unpolarized part 
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FIG. 5: Asymmetries for the difference of charged pion pro- 
duction with l^g{x) and —Ag(x). 



constraint power of experimental data are weaker than 
that of the tt'^ asymmetry below the region. However, 
the value of Rasym becomes larger than Rsta above the 
region. The asymmetry would have the same impact on 
Ag(x) as the 7r° production. Although more luminosity 
is needed in comparison with tt" production, it is useful 



TABLE I; The value of parameters q, -Rsta, and Rasy^ 



PT (GeV) 


9 


10 


11 


12 


13 


a 


0.82 


0.80 


0.78 


0.76 


0.74 


Rsta 


7.2 


6.4 


5.7 


5.2 


4.7 


Rasym 


5.0 


5.1 


5.3 


5.5 


5.6 



to determine effectively the behavior of Ag{x) with the 
sign. 



{a^^ — a'^ ) of the asymmetry in Eq. ijSJ). These contri- 
butions are not negligible. In particular, effect of the qq 
contribution in the polarized part appears as the differ- 
ence between the absolute values of asymmetries. Contri- 
butions of all sub-processes are taken into account, how- 
ever the q(f'^ — > q(f'^ and gg' — > q(( processes are negli- 
gible. The difference is due to the positive contribution 
of qq process. The asymmetry with — Ag(a;) is therefore 
suppressed. 

Next, we evaluate the experimental sensitivity of this 
spin asymmetry. We compare the statistical error of the 
asymmetry 5AJ^j^~^ with one of tt" production SAJ^j^. 
An expected statistical error SA^^j^~^ is given by 



IV. 



SUMMARY 



1 vT 



(12) 



where P is the beam polarization, a is the ratio of the 
number of event for tt" and 7r+: a — N'^ /N'^ . 
are obtained by the integrated luminosity C and the un- 
polarized total cross section : N'^ = £a^ . The 
ratio of these statistical errors can be obtained by 



Rsta 



1 1 



a 



(13) 



The parameter a has energy dependence, and decreases 
with pt increases. 

Table. represents the value of these parameters. 

Rasym = /^LL ratio of asymmctries. The 

7r+ — TT^ asymmetry is about 5 times larger than the tt'^ 
asymmetry. In the region < 11 GeV, the statistical er- 
ror becomes larger than the rate of the asymmetry. The 



In summary, we have investigated the correlation be- 
tween the behavior of the spin asymmetry for pion pro- 
duction and the functional form of Ag{x). The experi- 
mental data by the PHENIX indicates the negative asym- 
metry at low Pt, and motivate us to modify the func- 
tional form of Ag{x) drastically. In order to obtain neg- 
ative asymmetry, the functional form of Ag{x) requires 
some restraints. By modifying Ag{x), the slight nega- 
tive asymmetry can be obtained at low pT- Moreover, 
we have indicated the existence of the negative polariza- 
tion of Ag{x) which keeps the asymmetry to be negative 
in the whole px regions. The large negative asymmetry 
is inconsistent with the theoretical predictions by using 
Ag{x) from polarized DIS data. However, experimental 
uncertainties are large at present. It is premature to con- 
clude that the pQCD framework is not applicable to vr" 
production in polarized pp collisions. 

Uncertainty of the tt" asymmetry coming from the po- 
larized PDF's with DIS data is correspond to the current 
statistical errors by the PHENIX. These data have the 
same constraint power on Ag{x) as present DIS data. 
The future measurements will provide useful information 
for clarifying the gluon spin content. 

Furthermore, we have proposed the spin asymmetry 
defined by the difference of cross sections for 7r+ and tt~ 
production. We have discussed an impact of the asym- 
metry on determination of Ag{x). In the asymmetry 

, the gg processes are cancelled out, and qg pro- 
cess becomes dominant. Ambiguity of the fragmentation 
functions can be reduced. The behavior of the asymme- 
try is sensitive to the sign of Ag{x). One can obtain new 
probe for Ag{x) in pion production at RHIC. 
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